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A study of the photophysical properties of nitrogen-vacancy (NV) color centers in diamond na¬ 
nocrystals of size of 50 nm or below is carried out by means of second-order time-intensity photon 
correlation and cross-correlation measurements as a function of the excitation power for both pure 
charge states, neutral and negatively charged, as well as for the photochromic state, where the cen¬ 
ter switches between both states at any power. A dedicated three-level model implying a shelving 
level is developed to extract the relevant photophysical parameters coupling all three levels. Our 
analysis confirms the very existence of the shelving level for the neutral NV center. It is found that 
it plays a negligible role on the photophysics of this center, whereas it is responsible for an increasing 
photon bunching behavior of the negative NV center with increasing power. From the photophysical 
parameters, we infer a quantum efficiency for both centers, showing that it remains close to unity 
for the neutral center over the entire power range, whereas it drops with increasing power from near 
unity to approximately 0.5 for the negative center. The photophysics of the photochromic center 
reveals a rich phenomenology that is to a large extent dominated by that of the negative state, in 
agreement with the excess charge release of the negative center being much slower than the photon 
emission process. 

PACS numbers: 42.50.Ar, 42.50.Ct, 81.05.ug, 78.67.Bf 


I. INTRODUCTION 


With the development of photonic quantum cryptogra¬ 
phy and quantum information processes, there is need 
for reliable and easy-to-use single photon sources. Such 
sources have been developed in recent yearsi, like single 
molecules, colloidal or epitaxial semiconductor quantum 
dots, and color-centers in diamond^^— Here, we are in¬ 
terested in the latter, namely the NV center, formed by 
a substitutional nitrogen atom adjacent to a vacancy in 
the diamond lattice. NV centers have found numerous 
applications recently thanks to their unique physical pro¬ 
perties such as excellent photostability2i^“Z and long spin 
coherence times^ as well as to improved control over both 
their productioitS and physical initialization protocolsiiS 
NV centers can be made available both in ultra-pure 
bulk diamond^ and ultra-small crystals^ii Applications 
range from high-sensitivity high-resolution magnetome- 
tryi2.-i^, to fluorescence probing of biological processes 
2 ^ 1 ^, solid-state quantum information processing^l^, 
spin optomechanic a^^i^^ , quantum optics2^“— , nanopho- 
tonics^^Ti^l, and quantum plasmonics42r— 

NV centers can take two different charge states with dif¬ 
ferent spectral properties : the neutral center NV®, which 
has a zero-phonon line (ZPL) around 575 nm (2.16 eV), 
and the negatively charged center NV", which has a ZPL 
around 637 nm (1.95 eV)i^ In addition to ZPLs, the fluo¬ 
rescence spectra of both centers exhibit a broad and in¬ 
tense vibronic band at lower energy. A single NV center, 
as a single-photon emitter, is characterized by a second- 
order time-intensity correlation function that exhibits a 
photon antibunching dip at zero delay. In a first ap¬ 


proach, the photophysics of NV centers can be modeled 
by a two-level system with two photophysical parame¬ 
ters, the excitation rate, and the spontaneous emission 
rate. However, with increasing excitation power, the NV 
center, more particularly in the negative state NV", can 
experience distinctive photon bunching at finite coinci¬ 
dence time, in addition to the expected antibunching at 
zero delay. This can be accounted for within a three-level 
system with additional photophysical parameters to des¬ 
cribe photon decays to, or from, the additional shelving 
level. 

The aim of this paper is to give a detailed description 
of the intrinsic photophysics of single NV centers of both 
charge states in surface-purified^ nanodiamonds (NDs) of 
size around 50 nm, or below, as a function of the illumina¬ 
tion power from a GW laser. Understanding the intrinsic 
photophysics of NV centers is required before implemen¬ 
ting small fluorescent NDs in a complex electromagne¬ 
tic environment, such as practical single-photon devices, 
which will modify the photophysicsIt is also useful to 
the applications mentioned above as most of them exploit 
the single-photon emitter nature of isolated NVs. 

The statistics of both NV charge states has been studied 
previously in ND samples similar to those studied here 
and it was found to be size dependent, with a larger oc¬ 
currence (> 80%) of the NV" center over the entire size 
range from 20 to 80 nmi^ In the present study, we further 
observed that most of the single NV centers being detec¬ 
ted as neutral from their fluorescence spectrum at low 
illumination power (<0.5 mW) progressively gain with 
increasing power a photochromic character in that they 
also exhibit the NV" ZPL in addition to the NV’ ZPL. 
Only a few of NV^s remain purely neutral over the entire 
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power range, which we call pure NV^ behavior. On the 
other hand, NV centers detected as negatively charged at 
low power remain so with increasing power, which we call 
pure NV" behavior. We first focus our attention on such 
pure NV° or NV" centers. In addition to such behaviors, 
we found that some rare NVs can see their charge swit¬ 
ching between neutral and negative^^ already at the lo¬ 
west excitation power. We also describe the photophysics 
of such a photochromic NV over the same power range as 
the non-photochromic centers and show how both pho¬ 
tophysics can be linked. This allows us to find valuable 
information on the dynamics of photochromism. 

The paper is organized as follows. The experimental me¬ 
thods are described in Section II. The three-level model 
used to interpret the experiment is developed in detail in 
Section III. Section IV focuses on the experimental re¬ 
sults and on the extraction of the various photophysical 
parameters as a function of the excitation power for a 
NV” and a NV" center. Section V describes the photo¬ 
physics of photochromism in a single NV. A summary is 
given in Section VI. 



II. EXPERIMENTAL METHODS 

Preparation of the ND sample was achieved following 
a procedure reported previouslyi^i^ Commercial HPHT 
diamond nanocrystals are first irradiated using high- 
energy electrons, then annealed at 800 °C in vacuum to 
produce the fluorescent NV centers, and finally annealed 
in air at 550 °C to remove surface graphitic compounds. 
Colloidal dispersion in water and further sonication al¬ 
low us to obtain a uniform solution of NDs. We consider 
here NDs with a typical size of 25 nm or 50 nm deposited 
on a fused silica substrate to minimize spurious fluores¬ 
cence)^ Single NV centers were optically addressed using 
standard confocal microscopy at room temperature)^ A 
CW laser light (wavelength Aexc = 515 nm) falling wi¬ 
thin the absorption band of the NVs is used to excite 
the NV fluorescence. Excitation light was focused onto 
the sample using an oil immersion microscope objective 
with numerical aperture NA= 1.4. The NV fluorescence 
is collected through the same objective and is filtered 
from the remaining excitation, i.e., with wavelength be¬ 
low 532 nm, by a dichroic mirror and a high-pass filter. 
The collected fluorescence is subsequently sent either to 
a Hanbury-Brown and Twiss (HBT) intensity correlator 
(see below) or to a spectrometer. An example of fluores¬ 
cence spectrum is shown in Fig. [IJb) for the NV" case. 
The ZPL corresponds to the resonant decay (Fig. [IJa)) 
at A = 637 nm while the wide fluorescence side-band is 
the phonon replica. 

The second-order time-intensity correlation function 
contains the information on the classical versus quantum 
nature of light. It reads in the stationary regime 

( 1 ) 


Figure 1: (Color online) (a) Schematic of the two-level sys¬ 
tem explaining the emission spectrum of the NV center, g : 
ground state, e : excited state. The several thinner horizontal 
lines represent the coupling with the phonon bath in the dia¬ 
mond matrix. The green arrow represents the excitation and 
the red ones the optical decay at resonance (full arrow) and 
out of resonance (dotted arrow) (b) Typical emission spec¬ 
trum of an NV" center in a 25 nm ND. 


where P 2 {t + T\t) = P 2 (t| 0 ) is the conditional probability 
to detect a photon at time t + t knowing that another 
photon has been recorded at time t. This probability is 
normalized by the constant single photon detection rate 
Pi{t) = Pi{0)- For a classical source of light > 1, 

whereas the observation of an anti-bunching < 1 

is a clear signature of the quantum nature of light )^""— In 
particular, at zero delay (?^^^(0) = 0 for a single photon 
emitter—, which means that the probability to detect 
simultaneously two photons vanishes. 

In practice, an HBT correlator (Fig. [D(a))^i^ is used to 
measure g^^\ Here, the fluorescence of the NV center is 
sent on a beam splitter, which separates the signal in two 
equal parts sent to two avalanche photo-diodes (APDs) 
named APDl and APD2. The APDs are connected to a 
time-correlated single-photon counting module to build 
histograms of delays between photon events detected by 
the upper « Start » APDl and the lower « Stop » APD2. 
In order to avoid unwanted optical crosstalk between the 
APDs, a glass filter acting as a short-pass filter at 750 nm 
and a diaphragm are added in both branches/^ In the 
standard configuration, no bandpass filter is added to the 
setup, in contrast with the configurations used to study 
photochromism (Section V), so that it can be used for 
both charge states of the NV (provided that the NVs do 
not experience charge conversion, which is the case for 


g(2)(r)=P2(t + T|f)/Pi(f) 
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the selected NVs in the present work). 


their fluorescence spectra. 
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Figure 2: (Color online) (a) Schematic of the HBT cor¬ 
relator (b) Time-intensity second-order correlation function 
of a NV^ center excited at 200 /r W and presenting an 
antibunching at zero delay, (c) Time-intensity second-order 
correlation function of a NV" center excited at 5 mW 

and presenting a bunching at finite delay together with the 
antibunching at zero delay. SP 750 : short-pass filter at 750 nm 


A. Correlation function of a three-level system 


In order to give a quantitative description of the three- 
level model we first review briefly the interpretation of 
the function using Einstein’s rate equations. Starting 
from Glauber quantum measurement theory^ we have 




(: I{t + T)I{t) :) 

mtw 


( 3 ) 


where I(t) is the quantum operator equivalent to the elec¬ 
tromagnetic energy flow absorbed by a detector at time 
t and :: represents normal orderingj ^^i^^ Introducing the 
creation and annihilation photon operators (t) and a{t) 
in the Heisenberg representation, we obtain : 

f2i/_N {a^t)a'fit + T)a{t + T)a{t)) 

^ Kat(t)a(t))P ■ ^ ^ 


For a two-level system with ground state g and excited 
state e, the creation operator a^(t) at time t is to a good 
approximation proportional to the rising transition ope¬ 
rator \e,t){g,t\^, we deduce 


III. THREE-LEVEL SYSTEM : THEORETICAL 
MODEL 


(2). ^ ^ P{e,t + T\g,t) ^ p{e,T\g,0) 
^ ^ P{e,t) p(e,0) 


( 5 ) 


Fig. HDb) shows a typical g^^^ function measured for 
a single NV at low excitation power. The experimental 
data are compared with an equation stemming from a 
two-level model^i^^ : 

5(2)(r) = (2) 

where r and 7 are the excitation and spontaneous emis¬ 
sion rates, respectively. Within this model g^^^r) < 1, 
which means that the emitted light is non-classical at 
any delay r. However, at higher excitation rate this 
simple model (Eq. [5]) generally fails. This is particularly 
true for NV" centers. In this case, as shown in Fig. HDc), 
the g^^^ function includes a bunching > 1 feature 

at finite delays, superimposed to the antibunching curve. 
This kind of correlation profile, which contradicts Eq. [51 
calls for a third level that traps the electron, preventing 
subsequent emission of a photon for a certain time4^^— 
Therefore, for these delays, the correlation function is 
higher than one. Though well known for the NV" center, 
the existence of a shelving level is less documented for its 
neutral counterpart but has been invoked theoretically^ 
as well as in electron-paramagnetic resonance 4^ From 
now on, we will describe phenomenologically the color 
center dynamics by a three-level model and deduce the 
intrinsic photo-physical parameters of pure NV" and 
NV° centers (Section IV) and of a photochromic center 
(section V). These NV centers have been selected from 


where p(e, t-\-T\g, t) is the conditional probability for the 
NV to be in the excited state e at time t + t knowing 
that it was in the ground state g at time t. Like in Eq. [3] 
this probability is normalized to a single event probabi¬ 
lity, i.e., the probability p{e,t) for the NV to be in the 
excited state at the previous time t. The main interest 
of these equations is to link the probability of detection, 
as given by Eqs. 1 and |3l to the emission probabilities 
p{e,t + T\g,t) and p{e,t) defined by the rate equations. 
Therefore, can be completely determined if the tran¬ 
sient dynamics of the emitter is known. 

In the context of the NV center, we must use a two-level 
system with a third metastable state to explain the bun¬ 
ching observed in the correlation measurements!^ Since 
the previous calculations only considered the excited and 
ground states involved in the fluorescence process we will 
admit (see^^i for a discussion and justification) that 
these results still hold with a three-level system if we re¬ 
place the excited and ground states by the levels 1 and 2 
in the Jablonski diagram, respectively (see Fig. 3). Here 
we neglect the channels 1 to 3 and 3 to 2, because the 
system is not supposed to be excited at these transition 
energies, in contrast with previous models4^4^ There, 
channel 3 to 2 was taken into account and channel 3 to 
1 was neglected, because the quantum yield associated 
with the NV relaxation was supposed to be close to unity, 
while recent studies show that Q ~ 0.6 — 

Therefore, we here take into account four channels. This 
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leads to four unknown parameters : the excitation rate, 
ki 2 , the spontaneous emission, k 2 i, and the two parame¬ 
ters ^23 and ^31 of the additional decay paths involving 
the shelving level. ^21 being the only radiative channel, 
we can write the autocorrelation function as : 


where P 2 it) is the population of state 2 at time t. p 2 {oo) 
represents the asymptotic limit of P 2 {t) when the tran¬ 
sitory dynamics approaches the stationary regime. P 2 {t) 
can be obtained by solving the system of rate equations 
defining the three-level system presented in Fig. |3l With 



Figure 3: (Color online) Jablonski diagram of the three- 
level system, including the ground state (1), the excited state 
(2) and a metastable state (3). Only the allowed transition 
channels taken into account in the model are shown. 

Pi, i = {1,2,3}, the population of state i we can write 
the following set of equations : 

Pi = -ki2Pi -I- k2iP2 + fcaiPa, 

P2 = fcl2Pl - (^21 + fc23)P2, 

P3 = k23P2 - ksiPs, 

l=Pl-bp2+P3, (7) 

where pi means time derivative of Pi(t). The use of rate 
equations instead of Bloch equations is here fully justi¬ 
fied since in ambient conditions, the coherence between 
levels is decaying very fasti^ Eqs. 0 show that the sys¬ 
tem is necessarily in one of the states at any time. The 
steady-state analysis of these equations permits to find 
the explicit definition of the fluorescence rate R at which 
the system emits photons This definition allows us to 
explain the saturation behavior of the NV fluorescence, 
i.e., R tends towards a finite value for increasing excita¬ 
tion power. However, we are here looking for the time- 
dependent analysis to find p 2 (r) as a function of the kij 
coefficients. If we eliminate pa from Eq. [3 we obtain : 

Pi = -(fci2 -b k3i)pi + (fc21 - fc3i)p2 -b fcai, 

P2 = ki2Pl - (fel -b k23)p2- (8) 

The resolution of this pair of equations, with the initial 
conditions pi( 0 ) = 1 and P 2 ( 0 ) = 0 , leads to P 2 (r) = 
p( 2 ,r|l, 0 ) and therefore to the expression of 
Important approximations can be made to obtain a 


simple expression. More precisely, in addition to neglec¬ 
ting channels 1 to 3 and 3 to 2, we suppose that : 

{^ 21 , ^ 12 } {^23, fcsi}. (9) 

This is justified since, even if the third level is considered, 
the associated rates are supposedly very small compared 
to the singlet rates ki 2 and ^ 21 . We will see latter that 
this is not really true for the NV" center but that the re¬ 
sults obtained are actually robust and keep their meaning 
even outside of their a priori validity range. Within the 
above mentioned approximations the second-order corre¬ 
lation function reads (see Appendix A for mathematical 
details) : 

5 ( 2 ) (r) = 1 - /3e-T'i^ -b (/3 - l)e-T'=^ , (10) 

where the parameters 71 , 72 and /3 are defined through 
the relations : 


71 ~ fci 2 -b fel, 

( 11 ) 

, , ^12^23 

72 - K 31 + , , , , 

kl2 + k21 

( 12 ) 

/3^1+ 

k3i{ki2 + k2i) 

(13) 


B. Determination of the kij coefficients 

The aim of this sub-section is to determine the kij 
coefficients for the specific measured NV centers. The ht 
to the function allows us to determine 71 , 72 and 

/3. From Eqs. 11-13 we deduce : 


CM 

1 

II 

CM 

(14) 

r. J '2 

= J. 

(15) 

, 7i72(/3-1) 

^23 = 01, 

pkl2 

(16) 


We then have three equations for four unknown variables. 
A fourth equation is needed to solve the problem entirely. 
In the experiment, we also have access to the radiation or 
fluorescence rate R, measured in This rate is simply 
the average number of photons that the APDs collect per 
second. It represents, up to a multiplicative coefficient 
associated with the photon propagation in the setup, the 
probability for the system to be in level 2 , multiplied 
by the transition probability ^21 to relax (supposedly by 
optical means) to the ground state. We have : 

R = ^k2lP2{+00), (17) 

where ^ is the collection efficiency of the system once 
the NV center has emitted its fluorescence. Note that the 
same formula was used in refs 1 ^ 1 ^ with a different defini¬ 
tion of p 2 (+oo) because of a different Jablonski diagram. 
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Figure 4: (Color online) Time-intensity second-order correlation functions for a given NV° center (green lines) and 

NV" center (red lines). The excitation power is varied from 0.2 mW to 10 mW. The blue lines correspond to a three-level model 
fit. 


However, this has a physical meaning only if k 2 i is asso¬ 
ciated with a pure radiative decay. Whereas the model 
of refi ^^'^^ implied a unity quantum yield, the quantum 
yield in our approach is defined by 


Q 


^21 

fel + ^23 ’ 


(18) 


In our phenomenological approach the third level is thus 
assumed to absorb all of the non-radiative transitions 
letting fc 2 i be a pure radiative decay. Finally we remind 
that the probability p 2 {+oo) needed in Eqs. [T7] and [S] 
should be calculated in the asymptotic stationary regime, 
which can be obtained by canceling all pi in Eql8] : 

P2i+Oo) = - - - - - - --ry-- - - j- -r. (19) 

—/C21 + ^31 + (^21 + fc23)(l + ^ 31 /^ 12 ) 


Now four equations are at hand so that the system can be 
inverted to determine all coefficients. This calculation in¬ 
volves the numerical resolution through Cardano’s algo¬ 
rithm^ of a third-order polynomial as given in Appendix 
B. 


IV. EXPERIMENTAL RESULTS AND 
EXTRACTION OF THE fcp PARAMETERS 


We here consider two representative examples of NDs, 
one hosting a single pure NV’ center, and the second one 
hosting a single pure NV° center. The NV’-center ND is 
about 25 nm in diameter, whereas the one hosting the 
NV® center is about 50 nm in diameterThe func¬ 
tion was recorded for both NDs with different excitation 
powers Pexc in order to study the evolution of the 



Figure 5: (Color online) Evolution of the 71 , 72 and /3 pa¬ 
rameters as a function of the excitation power. Points are fits 
to the experimental functions using Eq. 10 and the lines 
are linear regressions. The green (red) points and lines corres¬ 
pond to the NV° (NV") center. Errors bars are shown where 
they exceed the symbol size. They represent the standard de¬ 
viation of each value calculated from the covariance matrix as 
given by the fitting routine. 


coefficients. 

Fig. [1] depicts the function for the two NV centers 
and for excitation powers Pexc ranging from 200 gW to 
10 mW. The experimental curves are fitted with Eq. 10 
taking into account the correction for the incoherent 
background light collected by the APDsi^i^i^^ The ex¬ 
perimental antibunching dip does not drop to zero due 
to this incoherent background, which modifies Eq. 10 as 
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5eS.('r) = + 1 - , where p = S/{S + B) 

contains the signal S and background B contributions 
from the NV fluorescence and the spurious incoherent 
light, respectively. By recording the average intensity 
from the APD directly on the NV and at a location close 
to it, we experimentally determine p and the fit parame¬ 
ters of Fig. |4] as explained in refsj^i^i^^. It is seen that 
the increase of Pexc only induces a very small bunching 
for the NV°, contrary to what is observed for the NV’. 
However, the anti-bunching dip narrows with increasing 
power in both cases. Our observations are collected in 
Fig. [51 which shows the power evolution of the fit para¬ 
meters 7 i, 72 and j3. 

The general trends seen in Fig. |3]are confirmed in Fig. [5] 
for both NVs since 71 , which is associated with the an¬ 
tibunching contribution, is clearly increasing with Pexc, 
a fact which is reminiscent from Eq. 2. Furthermore, it 
is seen that 72 is also increasing significantly for the NV’ 
which is clear signature of the third energy level. As far 
as the /3 parameter is concerned, it remains at a constant 
value ,5 ~ 1 for the NV“, while it increases up to /3 ~ 7 
for the NV’. This behavior agrees with its definition from 
Eq. 10. Therefore, when /3 is very close to 1, there is no 
significant bunching, and when /3 increases with power, 
the bunching turns on. 

Now that the three parameters have been found, the 
parameters can be traced back. This will be done in two 
steps. 


A. kij parameters with constant ^21 

In Eq. 17 the collection efficiency ^ of the optical se¬ 
tup must be known precisely to extract the various kij. 
However, since this can only be estimated, we will first 
calculate the kij parameters by assuming that fc 2 i does 
not change with the excitation power. This hypothesis 
is intuitive because the parameter k 2 i, i.e., the sponta¬ 
neous emission rate, is supposed to be solely governed by 
the Fermi’s golden rule, which in turn depends on the 
electromagnetic environment only. In order to determine 
the value of ^21 we observe that according to Eq. 14 we 
must have 71 = k 2 i at zero excitation since in this case 
ki 2 = 0. From the linear regression for 71 (Fig. [5]) we 
deduce 'y[^\Pexc = 0) = k^i = 0.052 ns~^ for the NV®, 

and 'y[~\Pexc = 0) = k^i = 0.046 ns~^ for the NV’ 
(here the (0) and (-) exponents refer to NV® and NV’, 
respectively). These constants give radiative lifetimes of 
T 2 i'^ = 19.2 ns and = 21.7 ns, which are consistent 
with previous reports (see for example refsi^*^). 

Thanks to Eqs. 14 to 16, we deduce the three other 
parameters as shown in Fig. | 6 | for both NVs. The first 
point to notice is that for both NVs, the ki 2 parame¬ 
ter increases linearly with the pumping rate from zero 
to a value exceeding fc 2 i for Pexc — 2 mW. Moreover, 
for NV°, we see that fcai increases, but keeps very small 
values compared to the other parameters, in agreement 


(a) 





Figure 6 : (Color online) Evolution of fcq parameters with 
power excitation when k 2 i is fixed for (a) the NV'^ and (b) for 
the NV’. Errors bars are shown where they exceed the symbol 
size. These error bars are estimated from those of the 71 , 72 
and P parameters, see Eqs. 14 to 16. 


with the assumptions made in Eq. 9. However, for NV’, 
the same parameters are no longer negligible compared 
to the set of ^ 21 , ^12 values. They even overtake k 2 i for 
Pexc > 6 mW. However, we emphasize that assuming ^21 
constant, if natural, is actually not fully demonstrated. 
In order to check how robust this hypothesis is, we will 
now try to approach the values and evolutions of the kij 
parameters that were obtained here by modulating the 
collection efficiency ^ using Eq. 17. 


B. Modulation of the collection efficiency p 

Now we use the fourth equation Eq. 17 to calculate 
the kij parameters. As already stated, we do not know 
exactly the p parameter but, as it turns out, slight varia¬ 
tions in p can produce significant changes in kij. To find 
the correct value of p, we adopt the following procedure. 
We let its value vary continuously and calculate the evo¬ 
lutions of the slope and Y-intercept of the linear regres¬ 
sions made with the obtained k 2 i{Pexc) traces. The slope 
should vanish because ^21 is assumed to be constant, 
whereas the Y-intercept should reach the value obtai¬ 
ned previously, i.e., the value of 71 at zero excitation. 
Therefore, we calculate the evolution of the kij parame¬ 
ters as a function of the collection efficiency p. The re¬ 
sults are shown in Eig. 0 where the constant horizontal 
curves depict the values to be reached by the slope and 
Y-intercept. Eor the NV'’ (Fig. |7| (a)), there is indeed a 
value of p where the two parameters reach the assumed 
values (blue vertical line). This gives = 0.77x 10“^, in 

agreement with a rough estimate of our setup collection 
efficiency taking into account the various optical compo- 
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nents. However, for the NV" (Fig. [3 (b)), it is seen that 
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(a) 



Figure 7: (Color online) Evolution of the slope and the Y- 
intersect of linear fit of k 2 i (Pexc) with the collection efficiency 
^ for the (a) and the NV" (b). Straight lines are the values 
obtained previously with k 2 i fixed. 


the symbol size in Fig. [5] (a), right panel), which confirms 
the very existence of the third metastable level for this 
charge state. Regarding the curves, it implies that 
the optical channel 2 to 1 is favored, which prevents any 
significant bunching. Furthermore, for the NV° the nar¬ 
rowing of the antibunching dip is due to the increasing 
excitation rate as for a two-level system, i.e., Eq. 2 (see 

refi2£) . 

The analysis of the kij parameters for the NV" is more 
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the parameters do not reach exactly the previous values. 
Yet, there is an optimum ^ for which the kij approach 
the previous values (blue vertical line). This corresponds 
to = 1.6 X 10“^, which differs only by a factor 2 from 
the case. This appears reasonable since the measu¬ 
rements were not carried out the same day (optical align¬ 
ments might be slightly different) and the collection effi¬ 
ciency ^ also depends on the unknown transition-dipole 
orientation in both NVs. 


C. Comparaison of the photophysics of NV centers 
in both charge states 

Fig. [5] depicts the kij{Pexc) curves deduced from the 
previous optimization. It is found that the order of ma¬ 
gnitude of the coefficients is the same as for imposed k 2 i. 
In particular, the relaxation rate fcsi is unchanged be¬ 
cause it only depends on the fit parameters, Eq. 15. 
For the two NV centers, the excitation rate vanishes in 
the absence of any excitation power and linearly increases 
with Pexc as it should. The main difference between both 
centers comes from the evolution of k^s, k^i and k 2 i- In¬ 
deed, for NV^, k 23 and ksi remain very small compared 
to fc 2 i, which is almost constant ( 1 /A: 3 i ~ 500 ns and 
1/^23 ~ 1000 ns if Pexc tends to zero). Therefore, the 
third level plays little role in the photodynamics of the 
NV° center. However, it is worth stressing that although 
very small, ^23 is not zero for NV° (errors bars are within 


Figure 8: (Color online) Evolution of the kij parameters 
with the excitation power taking into account the collection 
efficiency for the NV° (a) and the NV’ (b) centers. Errors 
bars estimated as in Fig. [6] are shown where they exceed the 
symbol size. 

involved. Indeed, ^23 increases very quickly to reach the 
order of magnitude of ^ 21 , while the latter is increasing 
as well (at zero excitation power we have 1 /A: 2 i — 24 ns 
and 1/^23 ~ 500 ns). The increase of ^ 23 , associated with 
non-radiative transitions, actually explains the growth of 
the bunching feature on the g^"^^ curves. Although the 
physical justification of this finding is beyond our pheno¬ 
menological treatment, it is likely that the variation of 
k 2 i and ^23 with excitation power is due to a change in 
the local energy environment of the NV center at high 
power, in particular because the efficient coupling with 
the phonon bath in the diamond matrix is expected to 
be temperature sensitive. Increasing the excitation power 
could thus correspond to an increasing effective tempe¬ 
rature, subsequently affecting the relaxation dynamics. 

It is worth pointing out that there is a limitation in the 
analysis done for the NV’ center. Indeed, the very fact 
that k 23 and {^ 21 , ^ 12 } reach the same order of magnitude 
contradicts the hypothesis made in Eqs. 9-13. Actually, as 
already mentioned, the results obtained are much more 
robust that could be anticipated at first sight. This can 
be figured out by relaxing the constraint of Eq. 9 as done 
in the detailed calculation presented in Appendix C. The 




























































Figure 9: (Color online) Evolution of kij parameters with 
power excitation taking into account the collection efficiency 
for the NV" in the approximation-free rate-equation model. 
Errors bars estimated as in Fig.|6]are shown where they exceed 
the symbol size. 


results obtained with the new rate-equation model are 
shown in Fig. IH] for the evolution of the kij coefficients 



Figure 10: (Color online) Evolution of the fluorescence rate 
(upper panel) and the quantum yield Q (lower panel) with 
the excitation power for the NV" and NV° center in the rate 
equation model free of any approximation. Errors bars for 
Q are shown where they exceed the symbol size. They are 
estimated from those of the kij parameters, see Eq. 18. 

of the NV" center. The obtained values are very similar 
to those of the approximate treatment, thereby justifying 
the previous results. For the NV° the coupling to the third 
level is very weak and the modifications (not shown) are 
even smaller. 

To complete the analysis we also computed the quantum 
yield evolution as given by Eq. 18 and compared with the 
evolution of the fluorescence rate. Fig. [10] confirms that 
within the Jablonski model sketched in Fig. 3 and in the 
considered excitation regime, the quantum yield of the 
NV“ center is approximately constant Q — 1, in agree¬ 
ment with the intuitive fact that the third level does not 
play a significant role in the dynamics. In contrast, the 
NV" quantum yield decreases dramatically with increa¬ 
sing excitation power from a starting (5~lto(5~0.5 


at high power. This entails the fact that the NV" dyna¬ 
mics is strongly dependent on the excitation power as 
discussed before. 


V. PHOTOCHROMISM 

Photochromism of NV centers has been reported 
in ensembles of NV centers in CVD diamond films 
under additional selective illumination^, with single NV 
centers in 90 nm NDs under femtosecond illumination, 
which results in the photo-ionization of the negative 
center to its neutral counterpart^, with ensembles 
of NV centers in type Ib bulk diamond at cryogenic 
temperatures under intense CW excitation^, and with 
a single center in natural type-IIa bulk diamond under 
CW illumination)^ In this last report, a special scheme 
of cross-correlation photon measurements was applied 
in the emission band of both charge states to show that 
the collected fluorescence in the NV° and NV" states 
were correlated and originated form a single NV defect. 
Several studies have reported that charge conversion 
within the NVs critically depends on the illumination 
conditions ,57,63,64 

A complete understanding of NV photochromism is 
lacking but a widespread view is that the optical excita¬ 
tion, either CW or transient, tunes the quasi Fermi level 
around a NV charge transition level, thereby inducing 
charge conversioni ^^i^^ This scenario has been reinforced 
recently by electrical manipulation of the charge state of 
NV ensembles and of single NVs by an electrolytic gate 
electrode used to tune the Fermi energy^ Depending on 
the sort of diamond studied, photochromism is thought 
to be favored by the presence of electron donor or 
acceptor defects, such as nitrogen, in the neighborhood 
of the NV center.— Recently, it was also found that re¬ 
sonant excitation of the NV*^ and NV" states in ultrapure 
synthetic Ila bulk diamond can induce reversible charge 
conversion in cryogenic conditions even at low power— 
This was taken as evidence that the charge conversion 
process is intrinsic in this sort of diamond, not assisted 
by an electron donor or acceptor state. The goal of 
this section is to give additional information on NV 
photochromism detected in surface-purified NDs, 25 nm 
in size, subjected to a CW non-resonant excitation of 
increasing power. By comparing the behavior of a single 
photochromic center to that of non-photochromic centers 
in the same illumination conditions as described above, 
we gain valuable information on the rich photophysics of 
photochromism. 

For the purpose of studying photochromism, we use 
two additional configurations of the HBT correlator 
that differ only by the set of bandpass filters added 
in the interferometer branches. These configurations, 
called NV"A, and NV®/" respectively, are shown sche¬ 
matically in Fig. [TTJ In contrast to Fig. Ha), these 
two configurations add selective bandpass filters in the 
interferometer branches. The NV"/° configuration in 
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Figure 11: (Color online) HBT configurations used to study 
photochromism. (a) is the configuration, where a band¬ 

pass filter (central wavelength 675 nm, bandwidth 67 nm) 
adapted to the negative NV" center is inserted in the upper 
start branch, whereas a bandpass filter adapted to the neu¬ 
tral NV’ (central wavelength 593 nm, bandwidth 46 nm) is 
inserted in the lower stop branch, (b) is the configura¬ 

tion where both bandpass filters are interchanged compared 
to (a). 

Fig. [11] (a) (respectively NV®/’ configuration in Fig. [11] 
(b)) uses a filter selective to the NV" (respectively NV”) 
fluorescence in the start (respectively stop) branch. 
Therefore, in the NV"/” configuration, a single photon 
emitted by a NV" center and detected in APDl gives 
the « Start » signal to the counting module, whereas 
a single NV” photon subsequently detected in APD2 
produces the « Stop » signal. The NV”/" configuration 
works in just the complementary way. Note that our 
NV"/” configuration is similar to the cross-correlation 
technique used in ref. — . These schemes turn out to 
be very powerful to study photochromism since cross 
correlation can be expected only if NV" and NV” 
photons originate from the same defect center. Note 
that related techniques have also successfully been 
applied to identify various excitonic species emitted by 
single semiconductor quantum dots, see for instancei^^— 


(a) (b) 



Figure 12: (Color online) (a) Fluorescence spectra of the 
photochromic ND nanoparticle taken at 3 different excitation 
powers of nominally 0.5 mW, 1 mW, 2 mW from bottom to 
top, respectively, (b) Antibunching curve at short time delays 
and an excitation power of 5 mW. 

With the above setup, we have located a few rare NDs 
hosting a single NV center that showed charge conver¬ 
sion at low excitation power already. In the following, we 
consider such a ND hosting a single photochromic NV 


Table I: Table summarizing the experimental parameters 

Pexc., R, ffexp.(O), S, B measured on the photochromic NV 

( 2 ) 

for the observation of the gfexp. function. The fit parameters 
71,72 and P using a 3 energy-level model are also given. 


center. 

The relevant spectra are shown in Fig. IT^ a) . It is found 
that both the NV” and NV" ZPLs are seen at any excita¬ 
tion power. The corresponding function measured in 
the standard configuration of the HBT (no bandpass fil¬ 
ter) is shown in Fig. lT^ b). It reveals a clear antibunching 
dip at zero delay. A precise analysis taking into account 
the background light (see Table 1) confirms this finding 
since at high excitation power we observe that the back¬ 
ground B increases significantly with respect with the 
NV fluorescence signal : p = S/{S + B) ~ 0.6 while at 
low power p ~ 0.9 — 1. With this value and using the 
formula gi%. {''')/ + 1 — 1 /p^ = 9^'^Kt) we deduce the 
actual value of 5*^^)(0) ~ q ,6 ,25,55,72 these results, 

a natural interpretation for the observation of the NV” 
and NV" ZPLs together with NV uniqueness is that this 
particular ND is subjected to photochromism. In addi¬ 
tion to the antibunching dip at zero delay, it is seen in 
Fig. [12] that ( 7 ^^^ exceeds 1 at longer delaysi^ In agree¬ 
ment with the previous sections, this is evidence for the 
presence of a trapping level and calls for a three-level 
description of the photochromic NV. The values of 71 , 
72 and /3 parameters used for fitting the gexp. function 
agree qualitatively well with those obtained for the NV" 
considered previously in agreement with the fact that the 
system is acting more like a three-level system. 

Although our study of NV photochromism is limited to a 
particular example we suggest that the dynamics of the 
system involves probably all energy levels of the NV" and 
NV” with some possible hybridization. It could be inter¬ 
esting to know, whether or not, the third level involved in 
the photochromic case is identical in nature to the third 
level of the NV". The role of the environment or of the 
radiation power— on the dynamics could be investigated 
in the future to clarify this point. 

Switching now to the cross-correlation regime we show in 
Fig.[T3]the result obtained for Pexc. = 3 mW (similar fea¬ 
tures, not shown here, were observed at other excitation 
powers). It is worth emphasizing that cross correlations, 
as described in e.g. refs i^^’^^'^^ , allow us to characterize 
the transitory dynamics between the two NV configura¬ 
tions. Using a formalism equivalent to the one leading 
to Eqs. 3-5 we indeed obtain in the NV"/” configuration 
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Figure 13: (Color online) Time-intensity second-order cor¬ 
relation functions measured for the photochromic NV center 
in two different HBT configurations, (a) is the cross correla¬ 
tion NV"/° configuration; (b) is the cross correlation 
configuration. The temporal window is [-100 ns; 500 ns] and 
the excitation power is 3 mW. 


sketched in Fig. [TT] (a) 

(2),. p(e,NV»,i + T|s,NV-,t) 

=- p(e.NV-,0 - 

where p{e, NV”, t -|- T\g, NV’, t) is the conditional proba¬ 
bility for the NV to be in the excited level e of the NV” 
state at time t + t knowing that it was in the ground 
energy level g of the NV’ charged state at the previous 
time t. This also means that a first photon emitted by 
the NV” was detected at time t while a second photon 
emitted by the NV’ is detected at time t -|- t . In a sym¬ 
metrical way using the NV”/’ configuration sketched in 
Fig. [TT] (b) we get 




p(e,NV’,t + r|g,NV”,t) 
p(e,NV”,t) 


( 21 ) 


with similar definitions as previously but with the role of 
NV” and NV’ inverted. 

The experimental results corresponding to these two 
configurations are shown in Figs. fTTT a) and (b), respec¬ 
tively. Here, it is also important to have r > 0 in the 
calculation leading to Eqs. |20| and in order to have 
a clear physical understanding. However, the electronic 
delay A = 100 ns included in the HBT correlator setup 
implies that sometimes even a photon emitted, say at 
time to) by a NV’ is recorded after a second photon emit¬ 
ted later by the NV” state, i.e., at to+r. This corresponds 
to the ’negative’ delay part of the graph, i.e.. Fig. |T3Ia), 
which is actually associated with the inverse dynamics 
NV”/’, i.e., Eq.jnH Here for clarity we did not subtract the 
delay from the abscises in Fig.[T31 Therefore for r > A in 
Fig. fT^ a) we have gi%.{''') = 5^/o(''’~^)P^ + l~P^ while 

we have g^eJp.ij) = g^^^_ (A - r)p^ -|- 1 - for 0 < r < A 
with p = S/{S B) as previously. In the same way for 
Fig. 113Kb) we have g^exp.ir) = gQ^/_(T - A)p'^ + I - p'^ 

for T > A and gi%XT) = - t)p^ + I - p^ for 

0 < T < A. 

As it is clear from the definitions these cross-correlations 


should present some symmetries. In the present case we 
used a three energy-level fit, i.e. Eq. 10, for the theoreti¬ 
cal functions, Eqs. [301 and HTj The parameters obtained 
to reproduce the data are the same for Fig. [T3| (a) and 
(b) up to an inversion between the ‘positive’ and ‘negati¬ 
ve’ delay for each graph. Therefore, we get for Fig. (Ha) 
7 i = 0.046 ns 72 = 0.01 ns ^ and /3 = 1.2 for r < A, 
while we have 71 = 0.052 ns~^, 72 = 0.01 ns~^ and 
/3 = 1.5 for r > A. For the second cross-correlation curve 
the parameters are identical but the roles of t > A and 
T < A are inverted as it should be. We observe that 
these values are very close to each other and also from 
the one obtained in Table 1 at the same excitation power 
Pexc. = 3 mW. This confirms that the system acts here 
mainly as a NV’ center. 

Interestingly, the last finding implies that after the emis¬ 
sion of a photon in the spectral fluorescence band of the 
NV’ the delayed emission of a second photon in the spec¬ 
tral band of the NV”, i.e., the conditional probability 
given by Eq. j^D] is also characterized by the dynamics 
of the NV’ contrary to the intuition. Such behavior was 
reported in refi^ for NV centers in bulk. In particular, 
in this paper the time dependence of the conversion NV” 
to NV’ process (and its inverse) was studied using pulse 
sequences. It was found that the relaxation from NV’ to 
NV” is a very slow process occurring with a decay time 
~ 1 ps. This agrees with our finding in Fig. [T3| since, 
even if the full dynamics of the photochromic NV cen¬ 
ter is expected to depend on the energy levels of both 
the NV’ and NV” centers, it is clearly the NV’ charac¬ 
teristics which dominate during the transition associated 
with Eq. |20| A similar qualitative analysis can be done 
in the NV” to NV’ center conversion. The small dissym¬ 
metry between the two parts of the curves for positive 
and negative delays results from the presence of a small 
NV” contribution to the dynamics during the NV” to NV’ 
transition which is absent in the NV’ to NV” conversion. 
Clearly, this complex charged/uncharged transition dy¬ 
namics would deserve systematic studies in the future. 


VI. SUMMARY 

To summarize, we have experimentally studied the 
fluorescence photodynamics of NV’ and NV” centers 
in diamond nanocrystals of 50 nm size or below using 
HBT photon-correlation measurements as a function 
of the excitation power. The dynamics was theoreti¬ 
cally modeled using Einstein’s rate equations and the 
transition probability rates entering the three-level 
model developed to analyze the data were deduced and 
used to infer a quantum efficiency to both charge states 
of the NV. It has been found that the shelving state, 
though present, plays a very small role on the neutral 
center in those small diamond crystals. The narrowing 
of the antibunching dip observed with increasing power 
for this center is a simple power effect that does not 
affect the near unity quantum efficiency. In contrast. 
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the negative center experiences a distinctive photon 
bunching behavior at finite delay that increases with 
increasing power. This reflects the increasing role of the 
shelving state for this center, which in turn diminishes 
the quantum efficiency from near unity at low power to 
approximately 0.3 at high power. We have also studied 
the dynamics of a photochromic center. It reveals a 
rich phenomenology that is essentially dominated by 
the negative face of the center. In the future, it would 
be interesting to determine how the presence of, e.g., 
a plasmonic structure affects the dynamics of NV 
centers^^ in any charge state, in particular if the role 
of the shelving state in the NV“ center as well as the 
neutral face of the photochromic center can be modified. 


VII. ACKNOWLEDGMENTS 


This work was supported by Ageiice Natioiiale de la 
Recherche, France, through the NAPHO (grant ANR-08- 
NANO-054-01), PLACORE (grant ANR-13-BS10-0007) 
and SINPHONIE (grant ANR-12-NANO-0019) projects. 


Annexe A 

Here we briefly summarize the derivation of the main 
equations used in sections III and IV. We start from Eq. 8 
written in a matrix form as 

= M-P + J, (Al) 

where 

a = -{ki2 + ki3 + fcsi), 
b = k2i- fcsi, 

C = ki2 — fc32, 

d = — (fc21 + k23 + ^ 32 ), 

and 



For the sake of generality we here keep all transition co¬ 
efficients kij allowed by the three-level model. Such an 
Eq.\M can be formally solved by defining a new vector 
q related to P by the equation P = • q. This leads 

to the new equation q = which has the general 

solution : 


q(t) = q(0) -I- f dre • J 
Jo 

= q(0) -f M-i • [I - • J (A3) 


where M ^ is the inverse matrix of M and where the 
initial condition corresponding to a system in the ground 

state at time t = 0 is q(0) = P(0) = f q )■ Therefore, 


we deduce : 

P(t) = ■ (q(0) 


M- 


[I - e-^“] • J) (A4) 


In order to give an explicit form to the solution we need 
to diagonalize the M matrix. Therefore, we write 


M n = 


Al 0 
0 A 2 


(A5) 


The eigenvalues Ai and A2 of M are easily obtained from 
the secular equation det[M. — AI] =0. They reads : 


Al — 


I 


A 2 = 


- - -y^[{a + dy -Met[M]] 


-A/[(a -I- d)2 - 4:det[M]\. 


From Eq. 


n” p = 


we easily obtain 

Alt Q 


0 




with 


and 


^ / I/A 


0 

V 0 1 /A 2 


n 


J 


(A6) 


(A7) 


(AS) 


(A9) 


In order to determine completely the solution we need to 


specify the transformation matrix 11 = 


Ui U2 
Vi V 2 


whose 


column vectors = 


with f = 1,2 are solutions 


of the eigenvalue problem M ■ = AiX^. These eingen- 

vectors are determined up to an arbitrary normalization 
and here we choose 


X,, = 


b 

Ai d 


(AIO) 


Using Eq. IA7I we therefore get : 

P2{t) = viae^^^ + U2£e'^^‘ -I- U 17 -I- U2(5. (All) 


We can further simplify the solution since from 
Eq. IA9I and the definition of M~^ we easily deduce 
a 
£ 


{V 1 V 2 ) 


= -(^'l^^2) 

V 26 ). After regrouping all terms we obtain 

P2{t) Via A,* , / Via 


, i.e., via + V 2 £ = —(uiy-l- 


P 2 (oo) 


Via + V 2 £ 


+ ( 


- I)e^=‘ -h 


Via -I- V 2 £ 


1 

(AI2) 


which up to notations is equivalent to Eq. 10 of Section 
III if we write Ai = — 7 i and /3 = —— 


vici 6 -\-V 2 e ’ 
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Annexe B 

The results obtained in Appendix A are exact and no 
approximation were made for calculating the coefficients 
7 i and /3. Now we will use the fact that {^ 21 ,^ 12 } ^ 
{^ 23 ,^ 31 } to simplify and explicit these coefficients. 
First, we point out that we have 

^ ^ (1 - Ji/li) ^ 7i72 
(71 - 72 ) Ji 

^ ( 1 -^ 31 / 71 ) det[M] 

+ dY — 4:det[M.]) ksi 

Therefore from Eas. lA6l and lBll we see that all coefficients 
7 i and /3 can be expressed as functions of a+d and (iet[M]. 
Up to the first order we have 


solution we use Eqs. 17 and 19 for the single-photon rate 
R. After insertion of Eqs. IB5I and IB7I into Eq. 17 and 
some lengthly rearrangements we finally obtain 

^12 “ 71^12 + 7‘(7i + T— )ki2 H—^ = 0 (B8) 

4 hi 4 

with X = {P — 1)7i 72//3- The three roots of this cubic 
equation can be obtained numerically using Cardano’s 
method. 

Annexe C 

We can generalize the analysis made in the previous 
Appendix B without any approximation. For this we first 
remark that we have the exact relations 


a -I- d ~ —(fci 2 -I- hi), 

det[M] ~ A: 3 i(fci 2 -b hi) + ^ 23 ^ 12 , (B2) 

Therefore, up to the same order, we have 
\/{{a + dY — 4det[M]) ~ —(a + d)[l —4det[M]/(a-|-d)^]. 
These lead to 


71 ~ -(a -b d) = fci 2 -b hi, 


72 - 


det[M] 


— hi + 


hak 


23«'12 


(B3) 


a -b d A:i 2 -b k 2 i 

which are Eqs. 11 and 12, respectively. Finally we have 


7172 = def[M], 

7 i +72 = -(a-bd), 

R ■ det[M] = 4^21^12^31 (Cl) 

which can be obtained after some manipulations from the 
definition of the M matrix and of Eqs. 17,19. If we now 
introduce the definitions 

A = det[M] = {ki 2 + k 2 i)k^i -b {ki 2 -b k 3 i)k 23 , 

B = —(a -b d) — k 3 i = ki 2 -b ^21 -b k23, (C2) 

C = R - det[M]/4 = ki 2 k 2 i 


— 1 det[M] 
(a -b d) k 3 i 


k23kl2 

k3l{ki2 -b /c 2 l) 


(B4) 


which is Eq. 13. 

In order to solve the system of equations 11-13 we first 
eliminate ^23 from Eq. 12 using Eq. 13, i.e. 


_ 7ifc3i(/3 - 1) 
ki 2 


(B5) 


Inserting this result into Eq. 11 leads to 

k3i = I (B6) 


which constitutes our first parameter solution. In order 
to determine the other parameters we insert the value 
obtained for fcai into Eg. IBSI and Eq. 11 to obtain 


we obtain after lengthly but straightforward calculations 

fci2 - Bki2 + C + A- Bk3i = 0. (C3) 

To complete the resolution of the system we use the exact 
relation Eq. IBII We finally obtain 

k 3 i = 727i/[7i - 72)18 -b 72], 
ki2 = B/2 - - 4{C + A - Bk3i)]/2, 

^21 = C/ki2, 
k23 = B — ki2 — k21- 

The minus sign was chosen in the second equation for ki 2 
(i.e. solution of Eq. IC3I) since a Taylor expansion at low 
excitation power when 72 ^ 71 , /3 ~ 1 gives for the two 
roots 


^23 = 


7172 (/3 - 1 ) 
fcl2/3 


k2i = 7i - ki2- 


(B7) 


k32 ^ 7 i/ 2[1 ± J(1 - I^)] 71/2 ± 71/2 T (C5) 


In other words, all parameters are now expressed as a In order to have the linear regime ki 2 oc i?/4 we must 

function of the excitation coefficient fci 2 . To complete the therefore impose the minus sign. 
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